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Abstract. Ultrathin polymer films on non-wettable substrates display dynamic features which have been 
attributed to either viscoelastic or slip efi'ects. Here we show that in the weak and strong slip regime efi'ects 
of viscoelastic relaxation are either absent or essentially indistinguishable from slip efi'ects. Strong-slip 
modifies the fastest unstable mode in a rupturing thin film, which questions the standard approach to 
reconstruct the effective interface potential from dewetting experiments. 

PACS. 83. 60. Be Linear viscoelatiscity - 47.50. +d Non-Newtonian fluid flows - 68.15. -l-e Liquid thin films 



Introduction. In recent years it has been shown that 
the physics of polymeric thin films on non-wettable sur- 
faces can be described, to an astonishing level of detail, by 
lubrication models derived fronr the Navier-Stokes equa- 
tion for simple Hquids [1,2]. However, ultrathin dewetting 
films exhibit unusual features in their rupture dynamics 
which show up in the morphology and velocities of dewet- 
ting holes [3]. It has been suggested that viscoelasticity 
plays an important role in these films, in particular when 
the polymer length scales become comparable to the film 
thickness. There is now a large number of modelling at- 
tempts to explain these features [4,5,6,7,8]. Most of them 
assume a generalized Maxwell- or Jeffreys-type dynamics 
for the stress-strain relation in these films, sometimes in 
combination with additional flow functions; all these as- 
sumptions are, while not entirely artificial, hard to soldily 
justify at present. In the absence of better knowledge, the 
Jeffreys-model therefore remains a useful starting point 
for modeling, with the idea to confront the predictions 
with experiment. 

As has been shown very recently, thin-film lubrication 
models can be classified into different slip-classes, and sep- 
arate models have to be derived for each class. There are 
models valid specifically in the limit of strong slip [9, 10] 
but also in intermediate slip regimes [10]. The distinction 
of different slip classes is essential for the description of 
dewetting experiments of PS-films on substrates with dif- 
ferent slip properties [11]. This last result has shown that 
slip effects can indeed explain the anomalies in the shape 
of dewetting films. 

As we demonstrate here, the distinction of different slip 
classes remains true for viscoelastic thin films of Jeffreys 
type. We show that it is easy to generalize the recently 



proposed lubrication model for Newtonian liquids in the 
strong-slip regime [9,10] to a Jeffreys model. We here ap- 
ply this model, as well as the recently developed model for 
the weak slip case [12], to determine the onset coirditions 
of rupture in unstable thin films. 

Model assumptions and lubrication equations. We be- 
gin with the bulk dynamic equations for the viscoelastic 
liquid. It is assumed as incompressible, hence the velocity 
field u — (ux,Uy,Uz) fulfills the mass conservation equa- 
tion 

V-u==0. (1) 

The equation of monrentunr conservation is given by 



Q 



d\i 
dt 



-VpR + V • r 



(2) 



where pr = p + V' is the augmented pressure, with p 
as capillary pressure and V' as the van der Waals-type 
dispersion forces. The traceless part of the stress tensor 
is described by a symmetric matrix r. Further, in eq.(2), 
d/dt is the total derivative, and V abbreviates the partial 
derivative vector with entries di, i = x,y, z. 

To complete the model we have to choose a consti- 
tutive relation for the stress tensor r. As argued in the 
introduction, we opt for the linear Jeffreys nrodel defined 

by 

{I + Xldt)T ^ Tj{l + X2dt)i (3) 

in which 7 is the strain rate, jij — diUj -\- djUi. The rates 
Ai, A2 govern the relaxation of the stress and strain rate, 
respectively. 

In order to derive the equations for a thin film of height 
z = h{x,y,t) we have, for the incompressible case, the 
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kinematic condition 



dth — — V|| • / dz u 



(4) 



where Vy = {dx,dy), and U|| = {ux,Uy). The boundary 
conditions at the free surface correspond to the vanishing 
of the stress tensor components tangential to the film sur- 
face (i.e., we neglect the vapor phase), while the normal 
component of the stress tensor obeys 



(r — prI) ■ n = 2crKn 



(5) 



where a is the surface tension of the film, 1 a 3 x 3 unit ma- 
trix, and K the local mean curvature with sign convention 
that K < for a spherical drop. Finally, in eq.(5), the nor- 
mal vector to the film is given by n = (— Vii/i, 1)/ ^/g, with 
g = 1 + (V||/i)^. The model is completed by the boundary 
conditions at the surface which are of Navier type, i.e.. 



u, = 0, 






(6) 



where b is the slip length. 

We now sketch the derivation of the lubrication model 
for strong and weak slip that can be derived from this bulk 
dynamics; the details can be found in rcfs. [12,10], and in 
the Appendix to this paper. 

First, we introduce a measure of the relative scale of 
the thin film height H, z = Hz*, to its lateral extension, 
L, {x,y) = L{x*,y*), and we define e = H/L < 1. Time 
is scaled by T = U/L, where U is the velocity scale. The 
stress tensor scales as 



— T- ■ 



(7) 



for {i,j) — {x,y) and, additionally, i ~ j. The remaining 
components scale as 



If' 



(8) 



The distinction between weak and strong slip lengths arises 
from the choice of balancing conditions between the forces 
acting on the film. In the weak slip limit, one has with the 
pressure scale P [12] 



while in the strong slip limit we need [10] 
PH 

^ £ . 



(9) 



(10) 



The Reynolds number Re = gUL/rj now scales as either 
Re = e^Re* (11) 



in the weak slip case, or as 



in the strong slip case, where Re* is the reduced Reynolds 
number of order one. Finally, the slip length b scales as 
b — 0(1) in the weak slip regime and b = (3se~'^ in the 
strong shp case. 

We first state the result for the strong slip case, details 
are given in the Appendix. Being interested here only in 
the conditions of thin film rupture, we restrict the discus- 
sion to the (laterally) one-dimensional case; the extension 
to the full two-dimensional case is straightforward. 

In the strong slip lubrication limit one ends up with 
the following system of equations (we put cr = 1), 



hRe*(dtu + ud^u) = hdx[dlh - V'{h)] + d^i^hq) - — 

Ps 



{l + \idt)q^{l + X2dt)dxU, 



dth + d^hu) =0. 



(13) 



where q is related to the stress tensor, see Appendix. Note 
that the system (13) readily reduces to the Newtonian case 
if Al = A2 = 0; the added complexity of the viscoelasticity 
is thus relatively minor in this limit. 

By contrast, in the weak-slip limit, one is able to derive 
the equation [12] 



(14) 



(1 + \2dt)dth + (A2 - Ai)9, \—Q^hR\ dth 



-dx 
where 

and 



h^ 



(1 + Xidt)'— + (1 + \2dt)bh' ) dxidlh - V'{h)) 



(l + \2dt)Q^-dx{dlh-V'{h)) (15) 



(l + \2dt)R = -hdx{dlh-V'{h)). (16) 



Re — eRe* 



(12) 



Note that for A2 -^ 0, eq.(14) collapses to a single equa- 
tion; this limit corresponds to the simplest Maxwell model. 
In the case Ai = A2 one recovers the thin-film equation of 
the Newtonian liquid with an extra multiplicative factor 
(1 + Xidt) on both sides. 

Linear stability analysis. We now turn to the linear 
stability analysis of a thin film which experiences a dis- 
persion forces which destabilizes it (i.e., V"{ho) < 0). The 
two different cases yield: 

A) Weak slip. The linear stability analysis is easily 
determined by assuming 

h = ho + Shi, Q^SQi, R = 5Ri, (17) 

where < (5 <C 1 with, in addition 

(/ii(x,i),Qi(x,t),i?i(x,t)) = (/ii,gi,i?i)e^'^"+'^*. (18) 

The resulting dispersion relation uj{k) can be expressed as 

{1 + \2Uj)lJ ==UJNil+Auj) (19) 
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where 



with 



^N{k)^-(^^+bhl^gik) (20) 



g{k)^k^ + eV"iho) (21) 

is the dispersion relation of the Newtonian liquid, and 



A = X2 + 



hi + mhl 



(22) 



From eq.(19) it is easy to sec that the structure of the 
dispersion relation of the Jeffreys film is identical to that 
of the Newtonian film. The range of unstable modes is the 
interval between the two zeroes of eq.(19) which is given 
by the two zeroes of ^(fc). Further, also the fastest unsta- 
ble mode is unaffected by viscoelastic relaxation. 

B) Strong slip. In complete analogy to case a) one puts 

h = ho + Shi , q = Sqi , u = 6ui (23) 

and, with eq.(18), one finds the dispersion relation 

(1 + XiUj){hoRe*uj + p-^)ui + 4/iofc2w(l + Aat^) 
+hlgik){l + Xiuj) ^ . (24) 

Again it is immediately evident that the range of unstable 
modes is unaffected by viscoelastic relaxation. The most 
unstable wavenumber km satisfies for the case Re* — 
which applies to the systems studied in [11] 

Ap.hlkt + hi {2kl + V"{ho)) I + ll^fj"^ = 0. (25) 

This result shows that the most unstable mode is strongly 
affected by slip, as was already observed for the case of 
a Newtonian liquid [9]. In addition we find that km also 
depends on the relaxation parameters Ai and A2. In the 
limit Al, A2 3> 1 eq. (25) simplifies to 



The derivation of the strong-slip lubrication model for 
the linear Jeffreys case follows closely both the calculation 
in the weak slip regime, and the strong-slip Newtonian 
case. As in [10], the starting point is the ansatz 

{u,w,h,pR,Tij) = {uo,wo,ho,pRo,Tija) (27) 

+ e^{ui,wi,hi,pjii,Tiji) 

where u and w are the velocity field components in x 
and 2;-directions, neglecting the transverse j/-direction. To 
leading order we find the equations 



TxzO = 

(1 + X2dt)d,uo - 
with the solution 

dzUa = c(x, z) exp(-i/A2) . 



(28) 



(29) 



(30) 



We select the solution c = since any other solution would 
correspond to a strong prestressing of the film at times 
t -^ —00. Therefore, uq = f{x,t), and from the mass 
conservation we have dxf ~ —dzW^, hence wq ~ —zdxf- 
It thus follows 



which reads in integrated form as 



(31) 



Tzzo = -Y- / dt'e^'-'">^^'{l + X2dt)dxf = -T,,o . 
-^1 J-00 

(32) 
To solve for f{x, t), we need to make use of the next order, 
which gives 

Re*{dt + fdx).f = dxTxxo + dzTxzi = -dxPm (33) 
where prq — ^dxxho — Tzzo ■ This can be written as 



± 



PI 
16 



V"{ho)p 



P 



Al 



PshoX2 



(26) 



Re*idt + fdx)f = dzTxzi + dxxxho 



(34) 



This result also holds for Re* 7^ 0, but the condition cor- 
responding to eq.(25) is much more involved. 

Conclusion. Based on the derivation of lubrication mod- 
els for thin-film dynamics of Jeffreys type we conclude that 
both in the weak and strong slip limits, linear viscoelastis- 
tic effects are essentially absent for film rupture. By con- 
trast, strong slip affects the the most preferred wavenum- 
ber, which now also depends on the relaxation parameters. 
In particular, from eq.(25) it appears that the standard 
approach for the reconstruction of the interface potential, 
which is based on the wavelength of the fastest growing 
mode is questionable for films subject to strong slip. 

Appendix A: Strong- slip lubrication limit for the Jef- 
freys model. 



+ —dx dt'e^'-''^/^'il + X2dt)dxf 

M j-00 

From the boundary condition at the free surface we find 
to second order 

{{Txxa - Tzzo) + Txzo{dxhQ)){dxha) = TxzI (35) 

and hence 

Txzi = -2(9j;/io)t2zo • (36) 

It remains to determine the second order result from bound- 
ary condition at the substrate. We have Txzi — f / P and 
can now integrate eq. (34) with respect to z across the film 
from to ft-o and obtain the system of eqs. given in the 
text (with / = u), and where 



Tzzn 



(37) 
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